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Abstract. In this paper, we analyze the quantum counting under the decoherence, 
which can find the number of solutions satisfying a given oracle. We investigate 
probability distributions related to the first order term of the error rate on the quantum 
counting with the depolarizing channel. We also implement two circuits for the 
quantum counting - the ascending- order circuit and the descending- order circuit - 
by reversing ordering of application of controlled-Grover operations. By theoretical 
and numerical calculations for probability distributions, we reveal the difference of 
probability distributions on two circuits in the presence of decoherence and show that 
the ascending-order circuit is more robust against the decoherence than the descending- 
order circuit. This property of the robustness is applicable to the phase estimation 
such as the factoring. 
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1. Introduction 

Since Grover demonstrated the database search algorithm [6], many applications of 
this algorithm have been studied [5, 7]. The quantum counting [4] is one of the most 
important applications with using the quantum Fourier transform. Given a quantum 
oracle, the quantum counting can find the number t of solutions of the oracle through 
in N elements with 0(\/tN) operations, whereas O(N) operations are required on a 
classical computer. The quantum counting is considered to be used for NP-complete 
problems because the quantum counting helps us determine the existence of a solution 
of these problems quadratically faster than a classical algorithm, solving whether the 
number of solutions is zero or non-zero. Recently, there have been works for numerical 
integrals based on the quantum counting, called quantum summation [1, 17, 8]. One 
can estimate the value of integrals with 0(l/e) operations for the desired accuracy e, 
while 0(1 /e 2 ) operations are needed by classical Monte Cario method. 

When we realize such quantum algorithms, decoherence is inevitable since our 
apparatus are surrounded with environment and open systems for us. Therefore treating 
the decoherence should be always one of significant problems. In the previous works, 
there have been only a few analyses of the decoherence, related to especially Shor's 
factoring [15] and Grover 's database search algorithm. Azuma [2] investigated the 
decoherence on Grover's algorithm by calculating quantum states in detail up to the fifth 
order term of the decoherence for a z errors. Shapira et al. [14] dealt with the decoherence 
on Grover's algorithm by changing the Hadamard transformation into some distorted 
operation. Yu et al. [18] represented disturbed unitary operations by the decoherence 
by focusing on quantum Hamiltonian and analyzed the evolution of quantum system 
on Grover's database search algorithm. Sun et al. [16] showed effects of environment 
based on the dynamic approach for quantum measurement through the example of 
the factoring. Several researchers have been investigated the decoherence by numerical 
calculations. Obenland et al. [12] simulated the circuits which factored the numbers 
15, 21, 35, and 57 as well as circuits that solved the database search for a trapped ion 
quantum computer. Niwa et al. implemented the general-purpose parallel simulator for 
quantum computing which could simulate not only the factoring and the database search 
algorithm [10] but also quantum error correcting codes [11], and revealed influences of 
the decoherence and another quantum error - operational error - to these algorithms. 
Although the decoherence on the factoring and the database search have been analyzed 
in a variety of ways, there have been no analysis of the decoherence on the quantum 
counting. 

In this paper, we investigate the decoherence related to the first order term of 
error rate on the quantum counting and expand our results to the phase estimation 
algorithm such as the factoring. The quantum counting is composed of two registers, 
called the first register and the second register. We assume the depolarizing channel as 
error models and calculate probability distributions on the quantum counting in two 
cases that the decoherence error occurs on each register. 
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We have another purpose in this paper to reveal which implementation for the 
quantum counting and one of the key quantum algorithms, the phase estimation, is 
robust against decoherence. Two quantum algorithms estimate a phase of a unitary 
operator by using the quantum Fourier transform and can be implemented in many 
ways by changing ordering of application of the unitary operations. In order to show 
the difference of effects of decoherence among different implementations, we implement 
two typical circuits - the as cending- order circuit and the des cending- order circuit and 
compare probability distributions on these circuits in the presence of decoherence. 

In the case that the depolarizing channel is applied on the first register in the 
quantum counting, we show that a probability distribution has many peaks caused by 
the decoherence at a distance of the power of two from correct peaks. We also show that 
probability distributions between the ascending-order and the descending-order circuit 
are the same. On the other hand, in the case of the second register, we reveal probability 
distributions on two circuits are completely different. In the ascending-order case, wrong 
outputs near correct one are obtained by the quantum counting, while in the descending- 
order case, two wrong outputs and N are obtained with high probability independently 
of an input oracle. Additionally, the correct output is obtained with higher probability 
on the ascending-order circuit than on the descending-order one. It follows that the 
ascending-order implementation for the quantum counting is more robust against the 
decoherence. 

We also show robustness against decoherence on the phase estimation with the 
ascending-order. Finally, we discuss weakness against decoherence on an efficient 
implementation for the phase estimation and the quantum counting, called the semi- 
classical implementation [13]. 

The rest of this paper is organized as follows: In Section 2, we begin by 
explaining the quantum counting and the depolarizing channel used in this paper, and 
implementing two quantum counting circuits. In Section 3, we analyze the decoherence 
on the quantum counting. We first consider analysis model for the decoherence, and 
then investigate probability distributions on the quantum counting in the presence of 
the decoherence on the first register and the second register. In Section 4, we also 
discuss influences of decoherence on the phase estimation and on the semi-classical 
implementation. In Section 5, we summarize our results. Finally, we give some detailed 
calculations of probability distributions in the case that the decoherence error occurs 
on the first register and the second one in Appendix A and Appendix B, respectively. 

2. Definitions, notations, and decoherence model 

First of all, we describe definitions and notations of the quantum counting used in this 
paper, and explain our decoherence model. We also implement two quantum counting 
circuits for revealing robustness against decoherence. 
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Figure 1. Circuit for the quantum counting 
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2.1. Quantum counting 

Suppose that there are unordered N := 2 n elements and a function / : {0, . . . , N — 1} — > 
{0,1} is given as an oracle (a black box). Let X be a set of elements which satisfy 
f(x) — 1 of N elements, that is, 

l (xel) 
(x(£l) ' 

and í := The goal of a quantum counting algorithm is to estimate the number t of 

solutions of the oracle / [4]. 

In order to estimate t, the quantum counting estimates the phase of a unitary 

operator G, called Grover operator, by using the quantum Fourier transform. Let 

\x) G Ti := (C 2 )® n , where (0 < x < 2 n — 1) be quantum states which represent 

2 n elements. Grover operator G is defined as G := U2U1 on 7i, where C/i := 

J2 x (-l) flx) \x)(x\,U 2 := 2\s)(s\ - l N ,\s) := ^EÍIo 1 !^- The operator G can be 

rewritten as a rotation on two-dimensional space. We divide the Hilbert space 7í into 

the "good" space H g :=span {\x)} and the "bad" space Hb :=span {\x)} and define two 

xel x<^x 
orthonormal states on each space: 



(1) 



... , 



a;SÏ 



The two-dimensional vector space spanned by the bases \b) and \g) is called Grover 
space. The Grover operator G can be represented on the Grover space as follows: 

_ . cos^ -sin^ . 

1 >\uü co S e ' ' 



where sin(0/2) := v^T^- 

Figure 1 shows a circuit for the quantum counting. We refer to the upper p qubits 
in this figure as a first register and the lower n qubits as a second register. The quantum 
counting algorithm is composed of the following five stages. 

(i) Prepare an initial state |o) 0( - p+n ^. 
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(ii) Apply the Hadamard transformation H to all qubits: ^= Xlm=o \ m ) ® \ s ) > wnere 
P := 2 P , \m) and |s) belong to the first and second register, respectively. 

(iii) Apply controlled-G to the second register according to the first register \m). 

(iv) Apply the inverse Fourier transform to the first register. 

p-i p-i 



m'=0 m=0 ^ ' 



m'=0 



(3) 



^ sin[7r(m / + /)] if sin[7r(m / - /)] 

Psin[7r(m' + /)/P] 1 ' Psin[7r(m' - /)/P] 1 

where / := P9/2n, |±) := ^(|&) =Fi |<7»- 
(v) Measure the first register and obtain rh for a good estimator of /. 

In order to determine the number of solutions from the measurement result rh, 
we calculate t := A^sin 2 (^/2) = iVsin 2 (7rm/P). Since the probability distribution of 
Equation (3) has peaks at rh ~ /, P — f(P 3> 1), we obtain the output t with high 
probability which is an approximation of t for the quantum counting. 

In this algorithm, we have to fix the parameter P that determines the precision. 
At first, by running this algorithm with setting P at v^ÏV, we obtain an approximated 
t such that |í — í| < 2n\/t + ir 2 . Then, by running it again with setting P at 20VtN, 
we obtain new estimation t as a more precise result. It is guaranteed to obtain t with 
probability at least 8/ir 2 [4]. It follows that the time needed for the quantum counting 
is 0(2?) = O(VtN)- 

In Figure 2, we show an example of the probability distribution of \m') in 
Equation (3) for the quantum counting, by setting p — 6, n = 8, X = {0, ...,12}. 
We also show a probability distribution of the corresponding output t' in Figure 3. The 
probability distribution of m' has two peaks, and that of t' has single peak at t that is 
the desired number of solutions. 



2.2. Decoherence model 

In order to analyze the decoherence, we need to model decoherence error. One of useful 
error models on classical information is the binary symmetric channel, which flips a bit 
with probability p and leaves it alone with probability 1 — p. We consider the following 
quantum decoherence model analogous to the binary symmetric channel, which is often 
used for analyzing error correcting codes. 

If we do not know anything about properties of errors that the quantum system 
suffers from, it is one of reasonable error models that quantum states are disturbed into 
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Figure 4. Two implementations of controlled-G™ operations. The ordcring (a) and 
(b) are called the ascending-order and the descending-order, respectively. Two circuits 
are equivalent if no error occurs. 



maximally mixed state as time goes on. In this paper, we assume that errors occur as 
local depolarizíng channel. 

Definition 1 (depolarizíng channel [9]). 

(1 -d)p + d- I - 

n j\ j a oPVo + <* x po x + o-ypo-y + a z pa z 
= {l-d)p + d , (4) 

where p is a density matrix on C 2 , a is the identity operator, and cr x ,a y ,a z are Pauli 
matrices. 

We apply this time-discretized error to each qubit at each unit time regardless of 
the existence of a quantum gate on the qubit. Since Equation (4) is represented by 
summation of classical events, we simulate this channel by applying a x , a y , a z , and a"o to 
each qubit with probability d/A] otherwise the state is left untouched. On our numerical 
calculations, we take an average of experiments by using only pure states. 



Theoretical Analyses of Quantum Counting against Decoherence Errors 



7 



2. 3. Two ímplementatíons for the quantum counting 

Generally speaking, unitary operations in circuit are not commutative. In the quantum 
counting circuit, however, all controlled-G gates are commutative. Figure 4 shows 
examples of equivalent implementations for the quantum counting, i.e. probability 
distributions on two implementations are completely the same in no decoherence case. 
One of dimculties in the analysis of decoherence is that we can not exchange the 
operators that are commutative in the system which is affected by decoherence. 

One of our aims in this paper is to reveal how different influences of the decoherence 
on equivalent circuits are and which implementation is more robust against the 
decoherence. In order to investigate it, we implement two typical circuits for the 
quantum counting. 

Definition 2 (two implementations for the quantum counting). We define two 
quantum counting circuits with the following ordering. 

ascending-order from the controlled-G 2 ' operation in Figure 4 (a). 
descending-order from the controlled-G 2P 1 operations in Figure 4 (b). 

The descending-order implementation is especially needed for an efficient 
implementation of the quantum counting and the phase estimation algorithm, which 
reduces the number of qubits on the first register to one. We discuss these efficient 
implementations in Section 4.2 later. 

3. Decoherence on the quantum counting 

In this section, we analyze influences of the decoherence on the quantum counting. We 
begin by explaining our analysis model and then investigate the decoherence on the first 
and the second register in two quantum counting circuits with the ascending-order and 
the descending-order. Finally, we discuss the robust implementation for the quantum 
counting against the decoherence. 

3.1. Analysis model for the decoherence 

The probability that the decoherence error occurs is considered to increase in proportion 
to the product of execution time and the number of qubits. As stated in Section 2.1, the 
quantum counting is performed in five stages: (i) preparing an initial state, (ii) applying 
the Hadamard transformation, (iii) applying controlled-G operations, (iv) applying the 
inverse quantum Fourier transform, and (v) measurements. The time needed for each 
stage is as follows: The stage (i) and (v) are considered to be one step. Application 
of the Hadamard transformations on the stage (ii) takes 0(1) and application of the 
inverse quantum Fourier transform on the stage (iv) takes at most 0(p 2 ). In contrast, 
application of controlled-G operations on the stage (iii) takes 2 P — 1, which is exponential 
to p. It follows that the decoherence error happens on the stage (iii) with exponentially 
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Figure 5. crj error occurs on the j-th qubit of the first register after application of 
controllcd-G* operations with the j-th control qubit. 



higher probability than on other stages. On our analyses, we restrict the position of the 
error only on the stage (iii). 

Suppose that the error rate d is small enough (d <C 1), then influences of the 
decoherence are approximated by the first order term of d. This means that the 
depolarizing channel is applied only once. Since the quantum counting has two registers, 
the register where the depolarizing channel is applied is either of the first register or the 
second one. 

Under these conditions, we calculate probability distributions on the quantum 
counting, and give some properties of the decoherence on two registers. Moreover, we 
compare influences of the decoherence on the ascending-order and the descending-order 
implementations. 

3.2. Decoherence on the first register 

We analyze the decoherence on the first register by considering the case that the error 
based on the depolarizing channel occurs once on the first register. Through the 
depolarizing channel, o~ x ,o~ y ,o~ z errors and identity operator a occur with the same 
rate from Equation (4). We first deal with the ascending-order circuit in Figure 4 (a). 

We investigate influences of the decoherence by calculating probability distributions 
on the quantum counting. For calculations, we need to represent the position of <7j error. 
As stated above, we already restrict the position of the decoherence on the state (iii). 
Let j and k be integers such that the error occurs 

• on the j-th qubit of the first register, 

• after application of controlled-G fc operations with the j-th control qubit, 

where < j < p — 1 , < k < 2^ . These parameters are sufncient for determining where 
and when <jj error occurs. The case that the error occurs on the j-th qubit before (after) 
controlled-G 2J can be represented by k = 0(k = 2 J ) respectively since our decoherence 
model is local. Figure 5 shows the position of o~i error on our analyses. The total number 
of applications of controlled-G operations before the error is 2 J ' + k — 1. 

We first focus on the position of peaks in a probability distribution on the quantum 
counting under the decoherence. If no error occurs, the probability distribution has only 
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two correct peaks near m ~ /, 2 P — f, as shown in Equation (3). Let Prob^ l ' : '' k \m') be 
the probability to obtain m' as a measurement result in the case of the above position 
of (Tj error. By calculation ProU l '^ k \m!) in Appendix A, we have 

Prob^ k \m') = 

i=0,x,y,z 

sin {n(m' + /) } 2 p ~ j sin {tt /2 p - j (mf + f) } 1 2 

sin {7r/2P-J- 1 (m' + /)} X 2p sin {7r/2P(m' + /)} 

sin7r(m'-/) 2 P ^' sin {n/2 p -i(m' - /)}! 2 

+ [2P-J- 1 sin {7r/2P-J- 1 (m / - /)} X 2^ sin {ir/2P(rrí - /)} J ' ^ ^ 

This equation has two strong peaks at m! ~ /, — / = 2 P — / that are the same as the 
correct peaks and weak peaks at a distance of ±2 p ~ j ~ 1 from the strong peaks, which is 
caused by errors. 

Proposition 1. The probability distribution related to the first order terra of the error 
rate on the quantum counting mainly has wrong peaks at a distance of the power of two 
from m ~ /, 2 P — f if the depolarizing channel is applied on the first register. 

We show the graph of ^Z^ i= Q iXi y^l^ P j^oProU t < j > k \m') in Figure 6 (a), which 
means the average of probability distributions in all error cases. We also show the 
graph of numerical calculations in Figure 6 (b). We did the experiments on Quantum 
Computation Simulation System (QCSS) [10] and took the average of 10 5 trials. On 
the numerical calculations, we set the error rate d = 4 x 10~ 3 so that the decoherence 
errors based on the depolarizing channel occur on the first register about twice on each 
trial. 

Proposition 1 follows the number of wrong peaks caused by the decoherence. 

Claim 2. The number of main wrong peaks is 0(p) in probability distribution related to 
the first order term of the error rate if the depolarizing channel is applied on the first 
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register. 

The probability distribution J2 i=Qxy z Prob^' k \m') does not depend on k that 
determines depth of <7j error, i.e. the time when the error occurs. 

Proposition 3. The probability distribution related to the first order terra of the error 
rate on the quantum counting is independent of depth of error if the depolarizing channel 
is applied on the first register. 

Proposition 3 means that there is no difference of influences of the decoherence 
between on the ascending-order circuit and on the descending-order one for the quantum 
counting, since controlled-Grover operators are commutative each other. 

Proposition 4. Probability distributions related to the first order terra of the error rate 
on the quantum counting are independent of the ordering of application of controlled-G 
operations if the depolarizing channel is applied on the first register. 



3.3. Decoherence on the second register 



We then deal with the case that the decoherence error occurs once on the second register 
in the quantum counting. In Section 3.2, we treat the error on the first register, which 
consists of control gates. In that case, error affects only the number of application m of 
controlled-G m . On the other hand, the error on the second register modifies the state 
on which Grover operator acts. Because the states \b) and \g) depend on the quantum 
oracle, we can not specify how the second register is disturbed by the decoherence. We 
need to consider the disturbance and action of Grover operator on a disturbed second 
register. 

We first show action of Grover operator on an arbitrary quantum state in order to 
deal with the operator on a disturbed second register. Any quantum state \<f>) G H is 
decomposed as follows by means of Gram-Schmidt orthogonalization: 

>W\b)+vW\g)+uí 



where u^\v^\u e r \v y e r ' G 
determined such that (b\e^) 
obtain the following lemma. 



(6) 



) G 7í b , \é?>) G Hg, and \e { f) and \ef) are 



14>)\ 



(g\ef ] ) 



0. By definition of Grover operator, we 



Lemma 5. For any quantum state \<p), Grover operator G can be rewritten as: 



G 



f 


cos 9 


— sin 9 





o\ 




sin9 


cos 9 
















-1 



















(7) 



satisfy |eH € H b , \e?>) G H g , (b\ef>) = (g\e 



3 W\ 

^9 



on four-dimensional space spanned by the basis states \b) , \g) , \e^), and \e g 9> ) , which 



0. 
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Figure 8. The probability distributions of the output t by numerical calculations 10 3 
trials with the error rate d = 4 x 10~ 3 . (a): on the ascending-order circuit, (b): on the 
descending-order circuit. 



Before disturbance, the second register is a superposition of the states \b) and \g). 
If some error occurs on the register, two states are disturbed into 

\b) ^u^\b)+v^\g)+uí^)+v?\eP), 

\g) -> \b) + \g) + u^le^) + v^\e^), (8) 

satisfying 

le^Ae^el·l,, \e^),\e^)en 9 , 

(b\eí ] ) = (b\eP) = {g\ef) = (g\e^) = 0. (9) 

Lemma 6. Any disturbed second register by the decoherence can be represented by the 
superpositions of \b) , \g) , |e^), |ejf^), \ef^), and \e^), satisfying Equation (9). Action 
of Grover operator is a rotation by 9 on two- dimensionat space spanned by \b) and \g), 
by 7r on |e^) and \e^), and by on \ef^) and |e^). 

Like the case of the first register, we consider that some error E, not necessarily the 
depolarizing channel, is applied on the second register after application of controlled- 
G k operations with the j-th control qubit, as shown in Figure 7. In this case, a 
probability distribution before measurement on the ascending-order circuit has peaks 
near m' ~ /, 2 P — / that are the same as ones in no error case, as detailed in Appendix 
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B. On the other hand, a probability distribution on the descending-order circuit has 
peaks not only near m! ~ f,2 p — f but also at m! = 0,2 p /2 with high probability, 
independently of the quantum oracle. By calculating an output t' = N sm 2 (nm' /2 P ) of 
the quantum counting, we obtain the following proposition. 

Proposition 7. The following wrong outputs of the quantum counting related to the 
first order term of the error rate are obtained with high probability if some decoherence 
error occurs on the second register: 

• Wrong outputs near t in the as cending- order case. 

• Wrong outputs and N in the descending-order case, independently of the quantum 
oracle. 

We show two graphs of outputs of the quantum counting with the ascending-order 
and the descending-order by numerical calculations in Figure 8 (a) and (b), respectively. 
These experiments were done 10 3 trials with the same conditions as the first register 
case except that the decoherence error occurs on the second register. 

The difference of positions of wrong peaks between two quantum counting circuits 
can be intuitively considered as follows: If a decoherence error occurs on the second 
register, influences of the error propagate to the first register by controlled-G operators. 
On the ascending-order circuit, Gs are applied from the controlled-G 2 ° with the Oth 
control qubit corresponding to the MSB of a measurement result. The influences 
therefore propagate the Oth control qubit (MSB) to (p — l)th control qubit (LSB). Since 
application of controlled-G m operations with the LSB needs more time exponentially 
than with the MSB, decoherence error occurs with exponential higher probability on 
controlled-G m operations with low control qubits. It follows that influences of the error 
propagate to only low control qubits. On the other hand, in the descending-order 
case, the influences propagate from the LSB to the MSB of the first register because 
of reversed ordering of application of controlled-G operations. Therefore not only low 
control qubits but also high qubits are affected by the decoherence. 

We finally consider the probability to obtain the correct output for the quantum 
counting. Peaks in the probability distribution on the descending-order circuit are 
distributed to four peaks whereas the probability distribution on the ascending-order 
circuit has only the correct two peaks. 

Proposition 8. The correct output of the quantum counting is obtained with higher 
probability on the ascending-order circuit than on the descending-order circuit if the 
decoherence error occurs on the second register. 

3.4- Robust implementation against the decoherence 

Proposition 4 shows that the probability to obtain the correct output on the ascending- 
order is the same as the probability on the descending-order one in the case of 
decoherence on the first register. Proposition 8 together with the proposition states the 
following robustness against the decoherence on two registers in the quantum counting. 
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Figure 9. Circuit for the phase estimation with the descending-order. Controlled 

with <j)j := e - 2m / v . 



uotations Rj are defined as Rj :=( * 



Claim 9. The ascending-order implementation for the quantum counting is raore robust 
against the decoherence than the descending-order implementation. 

4. Discussion on phase estimation algorithms 

In this section, we extend our results with respect to robust implementation against 
the decoherence to the phase estimation. We also discuss robustness of semi-classical 
implementation against the decoherence. 



4-1. Decoherence on phase estimation algorithms 

Phase estimation is one of key quantum algorithms, used in Shor's factoring [15]. 
Suppose that a unitary operator U has an eigenvector \u) with eigenvalue e 2nt(p , where 
the value of tp is unknown. The goal of the phase estimation is to find the phase 
ip. A circuit for the phase estimation is shown in Figure 9. The phase estimation is 
performed by application of controlled-?/ operations to the second register prepared 
to the corresponding eigenvector \u) initially. Like the quantum counting, the phase 
estimation has such equivalent implementations as the ascending-order circuit and the 
descending-order circuit, by changing ordering of application of controlled-?/ operations 
instead of controlled-G operations in the quantum counting. 

Claim 10. The phase estimation can find the desired phase on the ascending-order 
circuit with higher probability than on the descending-order one in the presence of 
decoherence on the second register. 

As mentioned in Subsection 3.3, influences of the decoherence error on the second 
register in the quantum counting propagate only to lower qubits of the first register with 
exponentially high probability on the ascending-order circuit, whereas these influences 
propagate to higher qubits of the first register on the descending-order circuit. This 
propagation of decoherence is applicable not only to the quantum counting but also to 
the phase estimation, though positions of wrong peaks caused by the decoherence are 
determined by a unitary operator U on the phase estimation. 
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Figure 10. Semi-classical circuit for the phase estimation. Rj is defincd with the 
results of previous measurements: R'- — , Q 
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4-2. Decoherence on efficíent ímplementations for the phase estimation 

An efïicient implementation for the phase estimation, called the semi-classical 
implementation, was demonstrated by Parker et al. [13]. We show the semi-classical 
circuit for the phase estimation in Figure 10, corresponding to the circuit in Figure 9. 
On this circuit, inverse quantum Fourier transform and measurements are applied as 
fast as possible so that only single qubit is required for the first register. By using 
semi-classical technique, Beauregard [3] reduced the number of qubits for the factoring 
about half as Ímplementations in Figure 9. It sounds so good for quantum computation 
since handling many qubits is considered to be difhcult. 

Here we discuss robustness of the semi-classical implementation against the 
decoherence. This implementation saves the number of qubits dramatically, although 
depth of circuits is almost the same as usual implementation. Since influences of the 
decoherence is considered to increase in proportion to product of the depth of circuits 
and the number of qubits, the semi-classical implementation is more robust against the 
decoherence from this point of view. 

We next focus on the ordering of application of controlled-?/ operations. On the 
semi-classical circuit, controlled-£7 operations must be applied from controlled-£/ 2ï> 1 
operations just as the descending-order case so that measurements are done as fast as 
possible. Since the semi-classical implementation has single qubit for the first register, 
most influences of the decoherence are caused by the second register. It means that the 
semi-classical circuit is less robust against the decoherence like the descending-order. 

We finally note the robustness of the semi-classical implementation for the quantum 
counting especially in order to solve NP-complete problems. Checking whether the 
number of solutions for these problems is zero or non-zero by the quantum counting 
helps us solve the problems. As we have shown in Proposition 7, the descending-order 
circuit has wrong peaks at or N independently of oracle. Since the semi-classical 
circuit is restricted to descending-order, it may not suitable for such problems. 

5. Concluding remarks 

In this paper, we focused on investigating influences of the decoherence related to the 
first order term of error rate on the quantum counting and revealing the difference of 
robustness against decoherence on two equivalent Ímplementations. 

In the analysis of decoherence on the first register, we showed that probability 
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distribution on the quantum counting had wrong peaks caused by the depolarizing 
channel at a distance of the power of two from correct peaks, and the probability 
distribution was independent of ordering of application of controlled-G operations. In 
the analysis on the second register, we first showed that wrong outputs were obtained 
near the correct one on the ascending-order circuit, whereas particular wrong outputs 
and N were obtained with high probability on the descending-order circuit. We 
then clarified that the correct output were obtained with higher probability on the 
ascending-order circuit than on the descending-order one. Consequently, the ascending- 
order implementation of the quantum counting was more robust against the decoherence. 

We also discussed the decoherence on the phase estimation. Similar to the quantum 
counting, the probability to estimate the desired phase by the phase estimation such 
as the factoring was higher on the ascending-order circuit. Moreover, we pointed out 
weakness of the semi-classical implementation against the decoherence because of the 
descending-order . 
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Appendix A. Probability distribution in the first register case 

We consider the case that the decoherence error based on the depolarizing channel 
occurs once on the first register. Suppose that the error occurs on the j-th qubit of the 
first register after application of controlled-G fc with the j-th. control qubit, as shown in 
Figure 5. Under the depolarizing channel, a x , a y , o z errors disturb each qubit with the 
same probability. We calculate probability distributions in each error case and take an 
average of probability distributions of all cases for simulating the depolarizing channel 
with pure states. 

Appendix A.l. cr x ,a y ,a z errors and identity a 

We first deal with a x error. Let ^2^Zo m j^ '■= in he indexes of the first register in the 
quantum counting. The quantum state before a x error is represented as follows: 



Since a x error flips a bit rrij of the quantum state \m), the state is disturbed by a x into: 




p-i 



1 P_1 

J2\m +(1-2171^) 

* m=0 




m 



mod 2?) 
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m=0 
P-l 



P-l 



F.T. 



m'=0 



m=0 l=± 



where |±) := ^(|6> ^i\g)) and ci ±} := ^ e 2ri ^e ±ií / 2 e ±i9 (" , -( 1 - 2 " 1 ^). 
By calculation of \cx^\, we have 



,(±)| 



1 P_1 

v m=0 



e 27riH^- e ±ie(m-(l-2m 3 )fc) 



sin[7r(m' ± /)] 



^2 sin[7r(m' ± f)/2P-^ 1 } 

2^'sin[7r(m / ±f)/2^\ 

X 



X cos 



n [ m'±f±2kíj ) /2"--> 



2Psin[7r(m / ±/)/2f] ' 
Similarly, we calculate \cf^\ where i — y,z, corresponding to <7j errors. 



,(±)| 



sin[7r(m' ± /)] 



| C W| = 



^2 2P-Í- 1 sin[7r(m' ± f)/Qp-i-*\ 

2 p -^' sin[7r(m' ±f)/2P-'] 
X 2Psin[7r(m'±/)/2P] 

1 sin[7r(m' ± /)] 

7f ' 2P-Í- 1 sin[7r(m' ± f)/2P-^ 1 ] 
2^sin[7r(m'±/)/2^] 



x sm 



7T ( rrí ± / ± 2k 



(A.l) 

/2P-J 
(A.2) 



x 



cos [n(rrí ± f)/2 p ~ ] ] 



x 



2Psin[7r(m'±/)/2P] 



(A.3) 



We also calculate \c^\ for symmetry of the depolarizing channel in Equation (4). 



,(±) 



1 sin[7r(m' + /)] 
Vl ' 2P sin[7r(m' + /)/2p] 
1 sin[7r(m' ± /)] 



^2 2P-J- 1 sin[7r(m' ± f)/Qp-i-*\ 
2P-J'sin[7r(m'±/)/2P-J] 



x 



sin [7r(m' ± /)/2 p ^] 



x 



2Psin[7r(m / ±/)/2P] ' 



(A.4) 



Appendíx A.2. Probabílity distribution of all errors 

Let Prob^'i'^ be the probabílity to observe m! as a measurement result, where i = 
x,y,z,0, < j < p — 1, and < k < 2 P — 1. Taking summation of probability 
distributions in all cases, we obtain the following probability distribution on the 
depolarizing channel: 

\ Y, Prob^\m>) = \ £ YM 2 



=0,x,j/,z 



=x,y,z,0 l=± 
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sin 7T (m 



' + /)] 
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2P-J-1 sin[7r(m' + /) /2P--?- 1 ] 
sin[7r(m' — /)] 



x 



2í , -^sin[7r(m / + /)/2^'] 
2f sin[7r(m' + f)/2P] 
2 p ^sin[7r(m' - f)/2 p - 



-i 2 



X 



1 V 

4 A^j=0,a;,: 



(A.5) 



2P-J- 1 sin[7r(m' - f)/2P~^- 1 ] 2Psin[7r(m / - f)/2P] 

z Prob^' k \m') does not depend on k. It means that this probability 
distribution is uniquely determined by the qubit where the error occurs, independently 
of when the error occurs. 



Appendix B. Probability distribution in the second register case 

Here we consider to calculate probability distributions of the quantum counting when the 
decoherence occurs on the second register. We deal with two probability distributions 
on the ascending-order circuit and the descending-order one. 



Appendix B.l. The ascending-order circuit 

We first consider the ascending-order circuit for the quantum counting. Suppose that 
the decoherence error occurs on the second register after application of controlled-G fe 
operations with the j-th control qubit, similar to the case of first register. 

Let r be the number of application of controlled-G operations according to the first 
register \m) before the decoherence occurs. 

r := krrij + (m mod 2 J ). 

The quantum state in the quantum counting before the decoherence is represented as 
follows: 

p-i 



Vp 



m=0 



cos ( r + - j 9 \b) + sin ( r + - ) 9 \g) 



By Lemma 6, the decoherence disturbs the state into 
p-i 

\m 
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?T5 S 
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where 
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satisfying Equation (9). By application of the rest of controlled-G operations and the 
QFT, 

p-i p-i 
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+e^(— ) („W| e W> + iu«| e ?>>) + (^)|ef ) + i„to>| e to>» }] .(B.2) 



We focus on peaks of each term in this equation. 



Appendix B.l.l. Peaks of the terms |±) 



Except global factors, a coeficient of the term |+) in Equation (B.2) is as follows: 

u (b) + ív (p) _ i^(ff) + iv (g)^ ^ 



2V2P 



p-i 
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M (6) + \V 
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2^2P 
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m=0 
P-1 



-e p > e ^ 

m=0 



Here we calculate only each summation instead of whole equation because we want to 
know where the term |+) have peaks. More precisely, we calculate the flrst summation 
p-i 

V.tti tnn I-77-I m' -4- t 1 1 

B.3) 



J^exp 

m=0 

and the second summation 

p-i 



27TÍ , , . 

— {m(m +/)} 



sin[7r(m' + /)] 



2Psin[7r(m / + /)/2P]' 



d 



(+) ._ 
ase •- 



^exp 



m=0 



2vri 



{mm' + (m — 2r)/} 



sin[7r(m' + /)] 



x 



2P-J-1 sin[7r(m' + Z)/^'" 1 ] 
2^'sin[7r(m / - f)/2?- j ] 



X cos 



^^+(l-|)/^/2- 



2Psin[7r(m' - f)/2P] ' 
The first summation has single peak at m! = —f = 2 P — f and the second summation 
has a strong peak at ml — f and weak peaks near m' — /, as shown in Figure Bl. It is 
easily seen that the term |— ) also has two peaks at m! = f, 2 P — f and weak peaks near 
m! = 2 P - /. 



Appendix B.1.2. Peaks of the terms \é^) and 
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We then consider peaks of the term |e^) which consists of two summations except global 
factors: 

p-i p-i 

5Z e ^ Í{mm ' +r/+(m ~ r) ^ } ' e^ {mm '~ r/+(m ~ r) ^ } , (B.5) 

m=0 m=0 

Since the term |e^) also have two same summations, we describe | e&) as terms |e^) 
and \ef S> ) for calculating peaks of the terms together. 
Let c?asc be the first summation in Equation (B.5). 
p-i 



"ase •- 



i ^ [27tí r , , , ,pi 

— >^ exp — - < mm + r / + [m — r) — > 

m=0 L k J 

sin[7r [m! + f )] 



X cos 



sin[7r (m' + f ) /2P-J- 1 ] 
2^sin[7rK + /)/2^] 
2Psin[7r(m' + /)/2P] ' 1 J 

Figure B2 shows a probability distribution of c4isc; which has a strong peak at m! = 2 p —f 
and weak peaks near the peak. The second summation in the term d^sc also has a strong 
peak at m' — f and weak peaks near the peak. 

Appendix B.1.3. Peaks of the terms |e^) and \e g g ^} 

The terms \ef^) and |e^), denoted by \e g ), have two summations: 
p-i p-i 

m=Q m=0 

By calculation of summations, we obtain two strong peaks at m' = /, 2 P — / and weak 
peaks corresponding to the peaks. 

Appendix B.1.4- Peaks of all terms 

In the case of the ascending-order circuit, all six terms in Equation (B.2) have only 
two strong peaks at m! = /, 2 P — f that are the same peaks in no error case and weak 
peaks near two peaks. Hence the overall probability distribution of the quantum state 
in Equation (B.2) is considered to have two peaks at m' = f, 2 P — f. 

Appendix B.2. The descending-order circuit 

We then investigate probability distributions in the descending-order case. Let r' be the 
number of controlled-G operations according to the first register |m) before the error. 

r' := m — 2 J mj — (m mod 2 JÍ ) + krrij. 

By simple calculation, we obtain the similar final state in Equation (B.2), which uses r' 
instead of r. We consider peaks of four terms \ef^), |e^), \ef^), and \e g 9 ^) in the equation 
because the terms |±) differ little from the case of the ascending-order. 



Theoretical Analyses of Quantum Counting against Decoherence Errors 
Appendíx B.2.1. Peaks of the terms \ef^) and |ejf^) 



20 



Like the ascending-order circuit, the term |e&) has two summations in this case by 
replacing r by r' in Equation (B.6). Let be the first summation of the term. 
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X cos 



7T l m! H r f 
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(B.8) 



^asc nas a strong peak at rrí = —P/2 = 2 p /2, as shown in Figure B3. This peak 
corresponds to —1 in the Grover operator in Equation (7). The second summation also 
has the same peak at 2 p /2. 

Appendix B.2.2. Peaks of the terms \ef^) and \e g 9 ^) 

Let c^gg be a summation of the term \e g ) corresponding to d^l in Equation (B.7). 
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(e g ) 
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sm 7T 



(m' + f )] 
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2 p ~i sm[7im'/2 p -:>] 



X cos 



tt ( rrí + -f ) IV 



2Psin[7rm72ï 5 ] ' v ' ; 

We show a probability distribution of d^^ in Figure B4, and obtain a peak at m' = 0. 
The second summation in \e g ) also has the same peak at rrí — 0. 



Appendix B.2.3. Peaks of all terms 



The overall probability distribution in the case of the descending-order has four strong 
peaks at rrí = f, 2 P — /, i.e. the correct peaks, and rrí = 0, 2 p /2, particular wrong peaks 
that do not appear in the ascending-order case and the first register case. The positions 
of these wrong peaks are fixed independently of an oracle and type of errors. 
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